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Abstract. It is known that, in continuous nicidia, composite solitons with hidden 
vorticity, which are built of two mutually symmetric vortical components whose total 
angular momentum is zero, may be stable while their counterparts with explicit 
vorticity and nonzero total angular momentum are unstable. In this work, we 
demonstrate that the opposite occurs in discrete media: hidden vortex states in 
relatively small ring chains become unstable with the increase of the total power, 
while explicit vortices are stable, provided that the corresponding scalar vortex state 
is also stable. There are also stable mixed states, in which the components are vortices 
with different topological charges. Additionally, degeneracies in families of composite 
vortex modes lead to the existence of long-lived breather states which can exhibit 
vortex charge flipping in one or both components. 



PACS numbers: 42.65Wi,42.82Et 
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1. Introduction 

Optical vortex solitons pLj, i.e., self-trapped beams containing phase singularities 
[21 El m [5l ini [7] , present an ideal setting for studying the relationship between topology 
and self-action effects, and may have applications to optical data transmission and 
processing [HIE]. However, in local nonlinear media vortex solitons are often destroyed by 
the modulational instability, and their orbital angular momentum (0AM) is transferred 
to muhiple filaments [Tnl[TniT2]. 

There are several approaches to suppressing the azimuthal instability of vortex 
solitons. Models with competing nonlinearities, such as cubic-quintic (CQ) [131 E], 
quadratic-cubic [151 [IS] nonlocal media [T71 [18] can support stable scalar vortex 
solitons. Two-component vortices may also be stabilized by the CQ nonlinearity [121 . An 
alternative approach is to apply spatial modulation to the self-defocusing nonlinearity, 
making its local strength to grow, with radius r, faster than [20|, or (more often) 
to apply a spatially periodic modulation to the refractive index of a nonlinear medium 
[21] . In the limit case of the deep periodic modulation, the medium reduces to an 
array of waveguides, the propagation through which is governed by discrete equations 
[221 123] . When the nonlinear self-action suppresses the discrete diffraction in the arrays, 
discrete solitons emerge [23]. The stability of discrete vortex solitons was predicted 
[211 [25] and subsequently observed in experiments [261 EZ] . Interestingly, the stability 
hierarchy is inverted with respect to continuous media: higher-order vortices, including 
"supervortices" (ring patterns built of compact discrete vortex solitons) ^2S] , tend to be 
stable, while their lower order counterparts suffer instabilities [23 [30l [311 [32] . 

On the other hand, vortex solitons can also be stabilized by the action of the 
cross-phase modulation (XPM) between two or more mutually incoherent components 
of a composite beam. In addition to the above-mentioned two-component model with 
the CQ nonlinearity [19], examples include multipole solitons [331 [31], two- and three- 
component necklace ring patterns [351 13SI [33 [3E], and their counterparts in the Gross- 
Pitaevskii equation for Bose- Einstein condensates [53 HD]- Interactions between two 
vortical composite solitons were studied in the model with the saturable nonlinearity 
|41j . Also studied were composite modes in which one component is vortical, while the 
other one is represented by a fundamental soliton [32j . In particular, solitons composed 
of vortex beams with oppositely rotating vortices, including symmetric counter-rotating 
pairs with the hidden vorticity (HV), whose 0AM is exactly zero, were predicted to be 
much more robust than their counterparts with the explicit vorticity [321 [33 [iQl [131 [11] ■ 
This feature can be demonstrated analytically in the framework of the one-dimensional 
two-component system, in which a counterpart of the HV states is represented by 
hidden-momentum counter-propagating wave pairs, with equal amplitudes and zero 
total momentum [H]. In the general case of an arbitrary number of symmetrically 
interacting components, the stability is determined by the total 0AM of the composite 
beam HH]. Vortex solitons of the HV type were recently observed in nematic liquid 
crystals [17]. 
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A specific example of composite {two-color) solitary vortices is provided by those 
in media with the quadratic {x^'^^) nonlinearity. While they are always unstable against 
splitting in the uniform media [101 [HI US], it was recently demonstrated that both 
single- and two-color vortices in x*^^) media can be stabilized by an external trapping 
potential [IH]. HV modes can be defined in terms of x^^^ systems too, assuming that 
the fundamental-frequency beam is built of two different components, corresponding 
to orthogonal polarizations, which are parametrically coupled to a single component of 
the second harmonic. The HV modes are unstable in that system (without an external 
potential), but the addition of the competing self-defocusing cubic nonlinearity makes 
them almost completely stable [IS]. 

This work aims to study combined effects of the above-mentioned approaches to the 
stabilization of vortices, viz., composite discrete vortex solitons in ring-shaped nonlinear 
lattices. We find that nonlinear modes of the HV type are subject to instabilities, which 
demonstrates that the above-mentioned inverse relation between the scalar and vortex 
stabilities and instabilities in discrete media, in comparison with continua, persists for 
composite modes as well. An additional feature of our discrete system is the existence 
of mixed-charge composite vortices. We conclude that their stability is tied to the 
topological charge of their brightest component. Further, an azimuthal modulation 
applied to discrete composite vortices may continuously deform them through a family of 
modes to discrete necklace solitons. When both components have equal total powers, this 
family is degenerate, and these deformations may be realized dynamically by perturbing 
a stationary mode. This can lead to simultaneous charge flipping of both components, 
similar to previous results in continua [1^. Additionally, due to the broken rotational 
symmetry of the discrete system, charge flipping of a single component only may occur 
too. 

We start by introducing the model, and obtaining analytically a class of "separable" 
nonlinear modes, in Section [2j Their stability is studied in Section [3j The dynamics 
of perturbed modes and vortex-charge flipping are presented in Section |4j Section |5] 
concludes the paper and discusses possible experimental realizations. 

2. The model and vortex modes 

We consider the discrete one-dimensional model governing the propagation of two 
incoherently coupled beams with amplitudes An{z) and Bn{z) through an array of 
nonlinear waveguides: 

id.An + An-l + A„+i + (|A„|2 + \Bn\^)An = 0, (U) 

id A + + + (|A„|2 + \Bn\')Bn = 0. (lb) 

We apply periodic boundary conditions, n-\- N = n, to define the ring-shaped arrays [cf . 
Ref . [Sni ED [S2] which considered in detail similar equations on an infinite chain] . The 
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system conserves the powers, Pa = J2n 1^™^, Pb = J2n l-^nP; the Hamiltonian, 



•(2) 



The current flows in the two components between adjacent sites, n and n + 1, are 
ji^n+i = 2Im(yl*yl„+i) and J^n+i = 2Im(i?*i?„_|_i). Discrete vortices correspond 
to circulation of the currents around the ring, and are characterized by the integer 
topological charges, 

'^A= arg(^n^n+i)' = 77- arg(5;5„+i), (3) 

n n 

which take values within the range of —N/2 < mA,B < N/2 [321 [5l]- There is no vortex 
in a component if its charge vanishes. 

Because the Hamiltonian ^ implies the Manakov-like nonlinearity, with equal 
coefficients of the XPM and SPM nonlinearities [53], it possesses an additional symmetry, 
associated with rotations that mix the two components, while preserving the power at 
each site: 

„ , , / cosv? sinv9 \ „ / cos6' isin6' \ , , 

\ — sm (f cos (fi I \ t sm tl cos u I 

where the rotation matrices act on vector {An,B„). The corresponding conserved 
quantity is commonly called the isospin, S = J^n^^B^, which exists in the continuum 
limit too [16]. Similar to the above definition of currents J^'t'+j, quantity S — S* = 2lmS 
is the total isospin-power fiow. Since the coupling between components is incoherent, 
i.e., they cannot exchange powers, this fiow must always be zero. Consequently, a global 
phase shift can always be applied to one of the components to set ImS* = 0. Therefore 
in this case the isospin has no physical significance, but the fact that it is conserved 
during the propagation will have some consequences later. 

We look for nonlinear modes as {An,Bn) = (f/„e*'''^^, V^e*^^^), where (3a,b are the 
propagation constants and f/„, Vn are the site amplitudes. Here we consider a special 
class of solutions, similar to necklace-ring vector solitons in bulk nonlinear media [35] , 
with constant total intensity on the ring, 

|f/nr + |Kp = / for all n. (5) 

Under this constraint, the nonlinearity ~ / is factorized, and the stationary equations 
for the amplitudes separate into two effectively decoupled linear equations, viz., 

(/ - f3A)Un + Un-l + Un+l = 0, (6a) 

(/ - + K-1 + K+1 = 0. (Qb) 

Exploiting the linearity of (6), we can apply the discrete Fourier transform, {Un, Vn) = 
^^(0^,63)6*®"'^, where = 2t[s/N is the phase winding of the Fourier mode with the 
discrete vortex of charge s. We thus obtain an analog of the dispersion relations [32] . 
/3a,b = 1 + 2 cos 9p^g, where the mode indices, p and g, may be different for two 
components with unequal propagation constants, Pa 7^ Pb- Note that they are 
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degenerate with respect to the sign of the mode indices p and q, thus we should take 
the superpositions, f/„ = a+e*®''" + a„e~*®''" and Vn = + 6_e~*®''", as a general 

solution. We can use the U{1) invariance of each component to set arga+ = — arga_, 
arg6+ = — arg6_ without the loss of generality. Applying the change of variables, 
a± = r-i-e^*^^, b± = s±e'^^^^, and A = (l/2)(/ — — — — s^), the constraint 
equation ([s]) takes the form of 

A = r+r_ cos[2(6pn + xi)] + cos[2(0<^n + X2)], (7) 

which must be satisfied for each value of n = 1,...,A^. This will fix some of the six 
parameters r±, s±,xi,2- Solutions with symmetries will have redundancies in the 
constraint equations, leaving free parameters. In the continuum limit, there is an infinite 
number of constraints, hence only symmetric solutions survive in that limit. On the 
other hand, additional modes may exist in the discrete system with a small number of 
sites. 

We will focus on a particular family of solutions with three parameters (P, 0, ip), 

Un = VP (cos 0) [(cos (^) e*®''" + (sin ^) e-^®^'"] , (8) 
Vn=VP (sin 0) [(cos 7) e^®'" - (sin 7) e"*®^'^] , (9) 

7 =^sin-i((cot2 0)sin(2^)), (10) 

which can be seen as a discrete generalization of the necklace-ring vector solitons [35] . 
Here P sets the total power of the mode, G [0,7r/2] determines the relative power 
in its two components, and (f G [—tt/Aj+tt/A] defines the azimuthal modulation. 
Notice that when the two components have different powers (0 7^ 7r/4), they require 
different modulation strengths, 7 7^ to maintain the constant intensity constraint 



(|5|. Obviously, (10) has a real solution only if \ip\ < (1/2) sin~^(tan^ 0). At the 
maximum allowed value of \(f\, the intensity of the second component is zero at some 
sites. By varying and if, a scalar vortex (0 = 0, = 0) can be continuously deformed 
into a symmetric composite vortex (0 = 7r/4, = 0) or a discrete necklace soliton 
(0 = vr/4, = 7r/4). Examples of these types of the modes are displayed in figure [l] 

In addition to solutions with the explicit {p = q) and hidden {p = —q) vorticities, 
for even A^ there also exist solutions with mixed vorticity, p = ±{N/2 — q). These appear 
because the constraint ([T]) only needs to be satisfied at a discrete set of points. 



3. The linear stability 

Linear stability is studied by introducing small perturbations of the form An = 

{Un + a^e^"^ + 6*e^*'^)e*^^'^, S„ = (K + c^e^'^ + (i*e^*^)e*'^s^ ^nd linearizing (1) to derive 
an eigenvalue problem for A. It is convenient to introduce vectors = {an,bn,Cn,dn), 
E„ = {An, A* , Bn, B^), writing the eigenvalue problem as 

{Hl + HNL)^n = iAdiag(l, -1, 1, -l)v„ (11) 
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Figure 1. Examples of different modes belonging to family (8)-(10) for the ring with 
iV = 6. (a) A scalar vortex {(j) = ip — 0) with charge p. (b) A symmetric vector vortex 
(0 = j,(p = 0) of charge {p,q). (c) A discrete necklace {(f> = j,(p — j). (d) Another 
discrete necklace, with = j,(p = —j. 



where Hi is diagonal in the component space (i.e., it does not couple a„ to 6„, etc. 
but couples adjacent sites: 

HLCln = (I — /^a)^?! + «n-l + 



and similarly for the other components (with (3a replaced hy (3b for c„, dn 
other hand, H^l couples different components, rather than different sites: 
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This can be compactly written as an outer product, H 



NL 



Et. 



Equations (11) can be solved numerically, as an eigenvalue problem of dimension 
AN . For large A^, it is relevant to consider some simple limits. For example, when ^9 = 
the problem simplifies significantly, the application of the discrete Fourier transform 
decoupling it into a set of "smaller" eigenvalue problems, each of dimension 4: 




where the matrices are defined as 



M 




(12) 
(13) 
(14) 
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Figure 2. Stable (shaded) and unstable (white) regions in the {P, (f>) parameter plane 
for composite vortices with different topological charges, indicated in the brackets 
{p,q), in the ring built of = 6 sites. Families of hidden- vorticity modes are 
represented by the top and middle panels in the right column. 



with Km,± = 2(cos9rrt±s — COS 9m), and |s| < N/2. In general, the eigenvalues cannot 
be obtained in an explicit form, but simple expressions are available in some limits. 
Stability diagrams of family (|8|-(10) in the limit of 99 = for = 6 are shown in figure 



[2] for different combinations of topological charges. 

When both components have the same charge (which corresponds to the explicit 
vorticity), the instability threshold is independent of 0, and the stability problem reduces 
to that for the scalar case: lower charges are unstable, while higher ones are stable [32] • 
Another situation takes place for the HV states with hidden vorticity: the HV mode 
(1,-1) encounters instability at a higher power than its scalar counterpart; ultimately, 
the (1,-1) mode becomes unstable at large powers, in contrast to stable HV solitons in 
continua [1^. Further, the HV configuration of the (2,-2) type also becomes unstable at 
large powers, in contrast to its explicit-vorticity (2,2) counterpart, which is completely 
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stable. 

Similar stability features are exhibited by the mixed vortex states of the (2,1), (2,- 
1) types: the state with the greater total charge p + q has a larger stability region. The 
stability of these solitons depends on the topological charge of the brighter component, 
i.e., whether the brighter component is stable in the scalar case. 

A qualitatively similar behavior is observed for other values of A^, with the stability 
dependent on whether the vortex charges are low (smaller than or high (larger 

than N/4). Vortices with the charges equal to N/A represent a special they 
form a one-parameter family of asymmetric vortex modes, introduced in [55j. This 
fact complicates the stability analysis, as the stability also depends on the additional 
parameter. 

When (j) = tt/A, an additional pair of zero eigenvalues appears for the hidden- and 
mixed vortex modes {p, —p), {p, N/2 — p), {p,p — N/2), and the linear stability analysis 
can no longer predict whether the modes are stable. The zero eigenvalues are associated 



with the degeneracy of family ([8j)-(10) with respect to (p. Calculating the values of 
the conserved quantities for the family, we obtain Pa = NP cos^ 0, Pb = NP sin^ 0, 
H = 2A^P(cosVcosep + sin^^cosGg) + A^PV2 and S = (l/2)A^Psin (20) [cos(v5 + 
7)5p,g + sin((/?-7)(5p,_g]. 

We see that for the HV mode with {p,—p), when = 7r/4, all quantities 
are independent of the azimuthal modulation which means that, under small 
perturbations, the input with (p = can cycle through solutions with different values of 
ip, hence we must consider the stability of the family as (p is varied. This family can be 
obtained by applying an isospin rotation to the ip = mode [46j : 



V I e-«^" ' ' ^^^^ 



which follows from the fact that (10) is solved by 7 = (y9 in this case. This rotation 



commutes with operator Hl + H^l in the linear stability problem, see (11). Following 
the result of |16], the stability is consequently independent of ip. This result is supported 
by direct numerical solutions of the linear stability problem. We stress that for other 
value of 0, the stability does depend on ip. 

The conserved quantities for the mixed vortex modes are independent of (p for any 
value of 0. However, the above reduction of the stability problem to the (p = case is 
not relevant, as the family cannot be generated by the isospin rotation. Therefore, the 
stability of the family depends on ip, as shown in figure |3] The azimuthal modulation 
lowers the instability threshold for the explicit vortex modes, and it can slightly increase 
the threshold for the HV mode. In practice, the family will become unstable when P 
exceeds the lowest instability threshold that occurs as ip is varied. 
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Figure 3. The stability of mixed vortex modes as a function of the azimuthal- 
modulation parameter (ip) and total power P. Shaded: marginal stability. Unshaded: 
linear instability 

4. Dynamics 

The existence of zero eigenvalues in the hnear stabihty problem means that higher-order 
terms will determine whether the degenerate families are stable. To check the stability, 
(1) was solved numerically using perturbed initial conditions. 

Figure |4] shows the propagation of a perturbed charge-2 HV mode in the waveguide 
ring built of six sites (A^ = 6). The high-frequency oscillations of site powers correspond 
to frequencies of stable eigenvalues, while the additional low-frequency oscillations 
correspond to the zero eigenvalue. During this oscillation, the power in each component 
acquires an azimuthal modulation, however the sum /„ = l^^p + l-B^P remains (on 
average) constant. 

When z ^ 25, the amplitudes of the m = 2 and m = —2 Fourier modes are 
equal and the topological charges of both components vanish [see figure |4]^b,d)]. Each 
component has two pairs of sites with equal powers, and the power vanishes at a pair of 
sites in one component. The beam profiles at this point resembles those in figure [T|c), 
which corresponds to a discrete necklace beam with (p = it /A. Increasing ip further, the 
charges of both components flip, and at 2; ~ 100 the azimuthal modulation vanishes, as 
if attains value 7r/2. There is a second charge flip at 2 ~ 150 {(p = Sir/ 4), this time with 
the other component hosting the pair of vanishing site powers. 

Thus we see that the slow oscillations correspond to an adiabatic cycling through 
the degenerate family of solutions parametrized by ip. No member of the family is 
subject to linear instability, therefore the oscillations persist indefinitely long (in excess 
of Az = 5000, according to our numerical results). 

In contrast to the picture for the HV model outlined above, instability can occur 
during the adiabatic cycling of the mixed vortex modes, since their stability depends 
on p. This is shown in figure [5} in which the initial condition is the perturbed p = 
mode. Initially, it experiences large oscillations in terms of its azimuthal modulation, 
while preserving its vortex charges and the power sums /„. At 2 ~ 400, an unstable 
value of p> is reached and an instability emerges, leading to irregular dynamics in which 
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Figure 4. Adiabatic dynamics of tlie degenerate iiidden-vortex family, with 
parameters P ^ = 7r/4, (p = 7r/8, topological charges (2,-2), and random 

perturbations added to the initial conditions at the 10% level. (a,c): On-site powers 
for the two components. (b,d): Vortex lines (red - positive charge, blue - negative 
charge) . 



neither the topological charges nor /„ are conserved. If, instead, the cp = tt/4 mode is 
used as the initial condition, this instability sets in immediately. 

An additional feature demonstrated by the discrete system is the broken 
conservation of the 0AM. Namely, the two components can exchange their angular 
momentum with the medium, as well as with each other. We show an example of this in 
figure |6| in which one component exchanges the angular momentum with the medium, 
leading to the periodic reversal of its topological charge, while the charge of the other 



component remains conserved. This type of dynamics is generated by solutions (|8j)-(10) 
with 7^ 7r/4 and chosen so as to make the value of 7 in (10) imaginary. The latter 



acts as a strong azimuthal perturbation, leading to the charge flipping of the second 
component. 



5. Conclusions 



We have studied the properties of composite vortex modes in circular arrays of nonlinear 
waveguides. The stability hierarchy of discrete composite vortices can be summarized 
by stating that the HV (hidden-vorticity) modes suffer instabilities above a critical 
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Figure 5. The adiabatic dynamics of the degenerate mixed- vortex family, with 
P — 1, (f) — 7r/4, (f = 0, topological charges (2,1) and random perturbations added to 
the initial conditions at the 5% level, (a) On-site powers of both components, (b) The 
total power at each site, /„ = -t- |-B„p. (c) The topological charges. 

power, while explicit vortex modes with high topological charges are stable. This 
hierarchy is opposite to that in continua. Additionally, mixed-vortex modes with 
different topological charges in the two components exist and can be stable. Degeneracies 
occur in these families of composite vortex modes, which results in long-lived breather 
states and persistent vortex-charge flipping. 

It should be stressed that the analysis was performed here for the relatively small 
ring chains, built of = 6 sites. For much larger rings, one may expect a change in 
the stability and dynamics, as, for a fixed diameter of the ring at — )■ oo, the system 
must go over into the continuum limit, with its reverse picture of the stability domains 
for the HV and explicit-vorticity modes. 

These effects are visible at high powers required for the self-localization in photonic 
lattices, therefore there is a possibility of observing them in experiments similar to those 
that were aimed at studying discrete vortices [261 EI] and multivortex solitons [56] in 
photorefractive crystals. A hexagonal lattice geometry corresponds to the ring built 
of 6 sites in our model. We can obtain appropriate experimental parameters for the 
observation of composite vortices from Ref. [31], which studied double-charge discrete 
vortex solitons, corresponding to the scalar (0 = 0) limit of our nonlinear modes. They 
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Figure 6. The charge- flipping dynamics of the single component, with P = 1/2. = 
7r/5, ip = 7r/10, topological charges (2,-1), and random perturbations added to the 
initial conditions at the 2% level. (a,c): On-site powers for the two components. (b,d): 
Vortex lines (red - positive charge, blue - negative charge). 



used a 20mm long crystal with a bias voltage of 2.2kV/cm, a lattice wave beam with 
a total power of 75fiW and a 30/im period. Linear propagation was observed with a 
probe beam at 532nm with a total power of P ~ 20nW, while the nonlinear regime was 
reached at P ~ 550nW. With these parameters, the propagation distance is long enough 
to observe the absence of discrete diffraction at high power (soliton formation) and the 
modulational instability of unstable modes. To observe composite vortex solitons, all 
that is required is to split the probe beam into two incoherent components, with the 
intensity and phase profiles generated using a spatial light modulator. Alternatively, our 
model can be realized directly in an integrated-optics setting, using a femtosecond- laser 
written ring of nonlinear waveguides and two incoherent beams [57], although because 
of the weaker nonlinearity, significantly higher beam powers would be required. 
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